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Abstract. A dynamical theory which accounts for all microscopic one-electron 
processes is used to study the spectral function of the ID Hubbard model for the 
whole (k, u>)-plane, beyond previous studies which focused on the weight distribution 
in the vicinity of the singular branch lines only. While our predictions agree 
with those of the latter studies concerning the tetracyanoquinodimethanc (TCNQ) 
related singular features in photoemission of the organic compound tetrathiafulvalenc- 
tetracyanoquinodimethane (TTF-TCNQ) metallic phase, the generalized theory also 
leads to quantitative agreement concerning the tetrathiafulvalene (TTF) related finite- 
energy spectral features, which are found to correspond to a value of the on-site 
repulsion U larger than for TCNQ. Our study reveals the microscopic mechanisms 
behind the unusual spectral features of TTF-TCNQ and provides a good overall 
description of those features for the whole (k, cj)-plane. 
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Most early studies of quasi- 

one-dimensional (ID) conductors such as tetrathiafulvalene-tetracyanoquinodimethane 
(TTF-TCNQ) have focused on the various low-energy phases, which are not metallic 
and correspond to broken-symmetry states [Tj. Recently, the resolution of photoemis- 
sion experiments has improved, and the normal state of these compounds was found to 
display exotic spectral properties [21 [3j 0] . However, such a metallic state refers to finite 
energies and thus is not described by the usual low-energy schemes. Based on a new 
quantum-object description of the ID Hubbard model [5], a preliminary version of the 
finite-energy dynamical theory considered here was used in the studies of Ref . [6] , which 
provided the weight distribution in the vicinity of well-defined branch lines only. Those 
lines refer to theoretical singular spectral features, which were found to describe the 
sharpest TCNQ related spectral dispersions observed in TTF-TCNQ by angle-resolved 
photoelectron spectroscopy (ARPES). Further details of such preliminary studies on the 
sharpest TCNQ related spectral features are provided in Ref. [TJ. In the mean while, 
a complete and more powerful version of the pseudofermion dynamical theory, which 
accounts for all one-electron microscopic processes, was presented in Refs. [HH]. (Such 
a general finite-energy theory recovers the known low-energy results in the limit of low 
energy, as confirmed in Ref. |10j.) 

In this paper the latter theory is suitably used to evaluate the momentum and 
energy dependence of the one-electron weight distribution of the model for the whole 
(k, a;)-plane. While both our predictions and those of Ref. jTT] (which refer to the 
TCNQ features only) agree with those of Ref. [6] for the TCNQ related singular features, 
we are also able to derive a theoretical weight distribution for the TTF related spectral 
features. For TTF the best quantitative agreement between the theory and experiments 
is reached for values of U larger than those preliminarily estimated in Ref. [6] for the 
TTF stack of molecules [12]. The U value found here for TTF is larger than for TCNQ, 
in agreement with results from the low-energy broken-symmetry phase pQ. Our study 
clarifies the microscopic processes behind the unusual spectral properties of TTF-TCNQ 
and provides a good overall description of its spectral features for the whole (k, a;)-plane. 
It also reveals that the electronic degrees of freedom of the normal state of quasi- ID 
metals reorganize for all energies in terms of charge and spin objects, whose scattering 
determines the unusual spectral properties. 

The structure of the quasi- ID conductor TTF-TCNQ consists of parallel linear 
stacks of planar TTF and TCNQ molecules [U El II]- Its partial charge transfer is 
0.59 electrons from the donor (TTF) to the acceptor (TCNQ) and thus the electronic 
densities are n = 1.41 and n = 0.59, respectively. Due to electronic correlations, the 
optical properties of the metallic phase depart significantly from Drude-free-electron 
behavior [1] and the finite-energy electronic structure, as probed by ARPES, deviates 
significantly from band-theory calculations [21 El H]. (See Fig. 7 of Ref. [1].) For 
energy values larger than the transfer-integrals for electron inter-chain hopping, the ID 
Hubbard model is expected to provide a good description of such correlations in quasi- 
1D conductors [U [13]. The model describes N = [N^ + iVjJ spin-projection a =|, j 
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Figure 1. Experimental peak dispersions (grey scale) obtained by ARPES on TTF- 
TCNQ along the easy-transport axis as given in Fig. 7 of Ref. [1] and matching 
theoretical branch and border lines. (The Z-point corresponds to the momentum 
k = 7T.) The corresponding detailed theoretical spectral- weight distributions over the 
whole (k, a->)-plane are plotted below in Fig. 2. While the theoretical charge-c" and 
spin-s" branch lines and c — s border line refer to the TTF spectral features found here 
(n = 1.41, U/t = 5.61), the charge-c, spin-s, and charge-c' branch lines correspond to 
the TCNQ dispersions (n = 0.59, U/t = 4.90) already studied in Ref. [BJ. 



electrons with densities n = N/N a and m = [AT-j- — iVj/iV in an ID lattice of N a sites. 
Except in the momentum axis of Fig. 1, we use units of lattice constant one, so that 
< n < 2. We denote the electronic charge by — e and define the Fermi momentum as 
kp = 7T 7T./2 for n < 1 (for electrons) and kp = 7r [2 — n]/2 for n > 1 (for holes). The 
model includes a first-neighbor transfer-integral t, for electron hopping along the chain, 
and an effective on-site Coulomb repulsion U. This is one of the few realistic models 
for which one can exactly calculate all the energy eigenstates and their energies [HIE]. 
Its low-energy spectrum belongs to the universal class of the Tomonaga and Luttinger 
liquid (TTL) theory [T5], [T6(, HQ]. However, its finite energy physics goes beyond (but 
is related to) the TLL and, until recently, it was not possible to extract from the exact 
solution the values of the matrix elements between the energy eigenstates. These are 
needed for the study of the finite-energy one-electron spectral- weight distributions. (The 
usual TLL theory can be considered as a special case of the finite-energy liquid used 
here, as we summarize below.) 

Our study focuses on the theoretical description of the unusual spectral features 
associated with the TTF stacks, which until now has remained an open problem. 
The model has both spin and i7-spin SU(2) symmetries. We denote by n and r] z = 
— [N a — N]/2 (and S and S z = —[Nf — -/VjJ/2) the 77-spin value and projection (spin 
value and projection), respectively. For U/t — > 00 all energy eigenstates correspond 
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Figure 2. Full theoretical distribution of the one-electron removal spectral- weight 
intensity (left) and corresponding line shapes (right). The figures include both the 
TTF related spectral features for n = 1.41; t — 0.35eV;[//i = 5.61 and those of 
TCNQ for n = 0.59; t = 0.40 cV; U/t = 4.90, respectively, as in figure 1. 

to electronic occupancies with fixed numbers of doubly occupied sites. However, the 
emergence of the exotic metallic state involves an electron - rotated-electron unitary 
transformation, such that rotated-electron double occupancy is a good quantum number 
for all U 1 1 values [5] . As the Fermi-liquid quasiparticles, such rotated electrons have the 
same charge and spin as the electrons, but refer to all energies and reorganize in terms 
of [iV a — N c ] 77-spin 1/2 holons, iV c spin 1/2 spinons, and N c spinless and ?7-spinless c 
pseudoparticles, where N c is the number of rotated-electron singly occupied sites [5]. 
We use the notation ±1/2 holons and ±1/2 spinons, which refers to the 7/-spin and spin 
projections, respectively. The ±1/2 holons of charge ±2e correspond to rotated-electron 
unoccupied (+) and doubly-occupied (— ) sites. The complex behavior occurs for the a- 
rotated electrons occupying singly occupied sites: their spin degrees of freedom originate 
chargeless a spinons, whereas their charge part gives rise to 77-spinless and spinless c 
pseudoparticles of charge — e. 

Based on symmetry considerations, we can classify the ±1/2 holons and ±1/2 
spinons into two classes: those which remain invariant under the electron - rotated- 
electron unitary transformation, and those which do not. The former are called 
independent ±1/2 holons and independent ±1/2 spinons, with numbers reading 
L c ,±i/2 = [v T Vz] an d L S} ±i/2 = [S =F S z ], respectively. The latter are part of r\- 
spin-zero 2z/-holon composite cv pseudoparticles and spin-zero 2z/-spinon composite 
sv pseudoparticles, respectively, where v = 1,2,... is the number of ±1/2 and 
— 1/2 holon or ±1/2 and —1/2 spinon pairs. The emergence of the exotic metallic 
state considered here involves a second unitary transformation, which maps the 
c pseudoparticles (and composite av pseudoparticles) onto c pseudofermions (and 
composite av pseudofermions) [E]. Such a transformation introduces shifts of order 
1/N a in the pseudoparticle discrete momentum values and leaves all other pseudoparticle 
properties invariant. As a result of such momentum shifts and in contrast to the c 
pseudoparticles and composite av pseudoparticles, the corresponding pseudofermions 
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have no residual-interaction energy terms [8]. 

As discussed below, the spectral-weight distribution of TTF-TCNQ is fully 
determined by the occupancy configurations of the c and si pseudofermions. We 
denote the latter by s pseudofermions. These objects carry momentum q. For the 
m = ground-state there is c pseudofermion occupancy for \q\ < 2k F and unoccupancy 
for 2k F < \q\ < n, whereas the s band is fully filled for \q\ < k F . (For m — 0, 
the exotic s band has a momentum width of 2k p.) The momentum, U/t, and 
n dependence of the a = c, s pseudofermion energy bands e a (q), group velocities 
v a{o) — de a (q)/dq, and Fermi-point velocities v c = v c {2k F ) and v s = v c {k F ) is 
provided by the exact solution j5j[H]. Under the ground-state - excited-state transitions, 
the ground-state a = c, s pseudofermions and holes are scattered by the a' — c, s 
pseudofermions and holes created in these transitions. Such zero-momentum-forward- 
scattering events lead to an overall phase shift Q a (q)/2 = Q° a /2 + {q'})/2, where 
Qt(q,W}) = 2E a '= c , s £ 3 ' a , (q,q>) AN a ,(q>), it <J> Q , a , (q, q') = -tt -q') 
is a two-pseudofermion phase shift whose q, q', n, and U/t dependence is provided by 
the exact solution, the momentum-distribution deviation AN a > (q') is that of the excited 
state, and Q° a /2 = 0,±tt/2 is a scattering-less phase shift whose value is well defined 
for each transition. The excited-state-dependent Fermi-point functionals Q c (±2kp) and 
Q s (±kp) fully control the spectral properties [8]. While the TLL theory involves a single 
interaction-dependent spectral parameter, 1 < £ < v^2 [H2, [H3, [10], the description 
of the finite-energy spectral properties requires interaction and momentum-dependent 
phase shifts n $ CjC (±2A;i?, q), n $ Ci s (±2fcp, q') n $ SjC (±fcF, q), and n $ SjS (±/ci?, q'), where 
2kp < \q\ < 7r for n > 1, < \q\ < 2k F for n < 1, and < \q'\ < k F . In the limit 
of low-energy, the scattering centers are created in the vicinity of the Fermi-points and 
our general description recovers the TLL theory with v c (±2k F ) = ±i> c , v s (±k F ) = ±t> s , 
and £o = 1 + [$ c , c (2fc F , 2k F ) - § c , c (2k F , -2k F )] = 2 [$ c , s (2k F , k F ) - <S> c , s (2k F , -k F )], as 
further discussed in Ref. [TO] . 

A crucial test for the suitability of the model is whether the observed ARPES 
peak dispersions correspond to the theoretically predicted sharpest spectral features. In 
figure 1 we plot the positions of the sharpest theoretical spectral features considered 
below but omit the corresponding detailed spectral-weight distribution over the (k, u>)- 
plane predicted by the theory, which is plotted in Fig. 2. The figure also displays the 
experimental dispersions in the electron removal spectrum of TTF-TCNQ as measured 
by ARPES in Ref. [3]. Figure 2 displays specifically the full theoretical distribution 
of the spectral-weight intensity (left) and the line shapes (right) corresponding to the 
same values of n and U/t as for the theoretical lines of Fig. 1. In the evaluation of the 
theoretical one-electron removal spectral features plotted in Figs. 1 and 2 we apply the 
improved pseudofermion dynamical theory of Ref. [8]. One of our goals is to find the 
value of U/t for which the theoretical weight distribution leads to the best agreement 
with the measured TTF related ARPES spectral features. 

The total number of ±1/2 holons (a = c) and ±1/2 spinons (a = s) reads 
M a> ±i/2 = L at ± 1 /2 + Y^Li u N au , where N av denotes the number of av pseudoparticles. 
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However, for the states which control the spectral properties of TTF-TCNQ, one has 
that N cv = for all v and N su = for v > 1 and for n < 1 (or n > 1), N c = N, 
L c ,+i/2 = [N a - N], and L c ,_ 1/2 = (or N c = [2N a - N], L c ,_ 1/2 = [N - N a ], and 
L c ,+x/2 = 0) and for m > (or m < 0), N sl = N t , L S;+ i/ 2 = [iV T - iVj, and £ s ,-i/ 2 = 
(or iV a i = iVf, L S) _i/2 = [iVj — iVj], and L S)+ i/ 2 = 0). Since the independent holons 
and spinons of these states are scatter- less objects [Sj, the weight distribution is fully 
determined by the occupancy configurations of the c and s = si pseudofermions. For the 
(k, u;)-plane regions of the ARPES data, the method used in our calculation involves 
specific processes associated with ground-state - excited-state transitions. For those 
which generate the dominant contributions one s hole is created and for densities n > 1 
and n < 1 a c pseudofermion and a c hole, respectively, is created. The low-energy 
TLL corresponds to processes where both such objects are created at momentum values 
in the vicinity of their Fermi points. Since the low-energy phase of TTF-TCNQ is not 
metallic and corresponds to broken-symmetry states, our results are to be applied for 
processes with energies larger than the gap, beyond the reach of TLL theory. 

For finite energy and U/t values all sharp spectral features are of power-law type, 
controlled by negative exponents. Important finite-energy processes are those where one 
a pseudofermion or hole is created at q away from the Fermi points and the second object 
is created at one of these points. The preliminary studies of Ref. [6] only considered 
such processes. They originate features centered on lines, u = uj a (q) = ±e a (q), in the 
(k, u;)-plane. In the vicinity and just below these lines, the spectral-function reads 012], 

B{k,Lu)^C a (q)(u a (q)-u)^. (1) 

When ( a (q) < 0, the spectral feature is a singular branch line. The exponent reads 
( a = —1 + (o(q), where Co(<?) is a functional whose values are fully controlled by the 
pseudofermion scattering. It reads, (o(q) = J2 L =±i^2a=c s^^ail) where 2A L a (q) = 
(lAN£ t + Qt( t( lFQ^ ( l)/^ 7! ') 2 ^ Qt( L( lFay ( l)/^ i s the scattering part of the overall phase 
shift defined above, lAN^ l = AqF a ^/[2Tr/N a ], and Aqp a ,i with t = ±1 is the Fermi 
point deviation relative to the ground-state values Lq° Fc = ilkp or iq% s = ik F . 
Expression ([1]) does not apply in the TLL regime, which corresponds to k « ±kp 
in Fig. 1, where the power-law exponent has a different expression [8]. The two regimes 
are separated by a small crossover region. Thus, the finite-energy normal state found 
here for TTF-TCNQ cannot be described by the usual TLL. 

There are other important types of finite-energy processes which were not 
considered in the preliminary studies of Ref. [6] . Those which involve creation of more 
than two quantum objects lead to very little weight. In turn, some of the processes 
where both a c pseudofermion or hole and a s hole are created at momentum values q 
and q', respectively, away from the Fermi points are important. When v c (q) 7^ v s (q'), 
such processes do not lead to singular spectral features and generate the background 
spectral weight all over the (k, cj)-plane, which although being in general small must be 
accounted for. (See Fig. 2.) Furthermore, a second type of sharp feature not considered 
in Ref. [6] corresponds to lines generated by such processes when both created objects 
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move with the same group velocity, v c (q) = v s (q'), and the spectral feature corresponds 
to a border line, u = uj BL {k) = [±e c (g) — e s (q')] o~v c (q),v a (q')i i n the (k, cu)-plane. The 
spectral function reads [9], 

B(k, LU) « C BL (k) (tU - LUBLik))- 1 / 2 , (2) 

in the vicinity and just above such a line. The TTF line called c — s in Fig. 1 and the 
weaker TCNQ bottom line of Fig. 2 are of this type. The latter weaker theoretical line 
is not plotted in Fig. 1, yet it clearly marks the lower limit of the experimental weight 
distribution. In the limit of U/t » 1, our one-electron weight distributions agree with 
those of Ref. [Hj. 

By careful analysis of the k, lu, and U/t dependence of the obtained theoretical 
weight distribution, we find that for n — 1.41 the electron removal spectrum calculated 
for t = 0.35 eV and U = 1.96 eV (U/t = 5.61) yields the best agreement with the 
TTF experimental dispersions. (The U/t = 5.61 TTF value is much larger than that 
preliminarily estimated in Ref. [6].) Remarkably, the only fitting parameter is U/t. For 
the considered values of n and U/t, the singular charge-c" and spin-s" branch lines and 
the singular c — s border line of Fig. 1 correspond to the sharp spectral features of 
the model one-electron removal spectral function for the ARPES (k, a;)-plane region. 
The fading parts of the theoretical charge-a = c', c" branch lines of Fig. 1, not seen in 
the experiment, correspond to values of the momentum where the constant C a (q) of the 
expression (CQ) in the vicinity of these lines is small. Although there is a reasonably good 
overall quantitative agreement between the theoretical one-electron weight distributions 
and the TTF related features measured by ARPES, there are apparent differences in 
the finest details, for instance the broadening of some of the sharp features predicted 
by the theory. However and in spite of the recent improvements in the resolution of 
photoemission experiments [2} [3], H], it is difficult to measure the weight-distribution 
finest details experimentally, in part due to the extrinsic losses that occur on anisotropic 
conducting solids [H]. Hence while our theoretical description provides the dominant 
microscopic processes behind the overall unusual spectral-weight distribution observed 
in the real material, it is difficult to judge which other smaller effects may play some 
role in the weight-distribution finest details. 

Our general study refers to the whole (k, u;)-plane and confirms the validity of the 
predictions of Ref. [6] for the TCNQ related spectral features: for n = 0.59 the finite- 
energy-electron-removal spectrum calculated for t = 0.40 eV and U = 1.96 eV (U/t = 
4.90) yields an almost perfect agreement with the TCNQ experimental dispersions, 
which correspond to the spin-s, charge-c, and charge-c' branch lines of Fig. 1 and the 
weaker border line shown in Fig. 2. There we plot the full theoretical distribution of 
the weight intensity resulting from electron removal both for n = 1.41; U/t = 5.61 and 
n = 0.59; U/t = 4.90 and the corresponding line shapes, respectively. 

Our results reveal that the ARPES peaks refer to separate spin (s hole) and charge 
(c pseudofermion or hole) objects for the whole energy bandwidth, whose line-shape 
depends on the interaction. For the Hubbard model, such a spin-charge separation 
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persists in the limit of low energy, where the quantum liquid becomes a TLL. An 
important exception is the TTF singular border line named c — s in Fig. 1 (and the 
weaker TCNQ border line shown in Fig. 2), which refers to a charge and a spin object 
moving with the same velocity. Before merging this line, the charge-c" branch line refers 
to q values such that |v c (9)| > v s , whilst for q and q' obeying the relation v c (q) = v s (q f ) 
such that < |f c (<z)| = l^sl*?')! < v s the singular c — s border line emerges. Such a 
feature does not exist at low energy because of differing charge and spin velocities. 

The transfer-integral values obtained here are about twice as large as those found 
by band theory, consistently with the experimental bandwidth being much larger than 
predicted by traditional estimates [2]. Moreover, our values for U/4t are of the order of 
unity and larger for TTF (U/t = 5.61) than for TCNQ (U/t = 4.90), consistently with 
the TTF-TCNQ broken-symmetry-states and optical properties [lj . The effects of the 
onsite repulsion U considered here for the electrons inside the solid, lead to a weight 
distribution that agrees with the general ARPES spectrum structure over a wide range 
of finite energies. (This is in contrast to band-theory calculations, as confirmed in Fig. 
7 of Ref. [1].) However, when the photoelectron is in the vacuum above the crystal it 
may create excitations in the substrate via long-ranged interactions beyond our model. 
The resulting inelastic losses as well other effects of finite temperature and long-range 
Coulomb interactions [19] are expected to be the mechanisms behind the broadening of 
the singular features predicted here [18]. Such effects lead to the broad peaks observed 
in the ARPES of Ref. [4] but do not change their overall distribution over the (k, u) 
plane, which remains as found in this paper. 
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